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ABSTRACT 

We use new interior models of cold planets to investigate the mass-radius relationships of solid exoplan- 
ets, considering planets made primarily of iron, silicates, water, and carbon compounds. We find that 
the mass-radius relationships for cold terrestrial-mass planets of all compositions we considered follow a 
generic functional form that is not a simple power law: log 10 R s = &i + 1/3 log 10 (M s ) — k 2 M^ 3 for up to 
M p « 2OM0, where M s and R s are scaled mass and radius values. This functional form arises because 
the common building blocks of solid planets all have equations of state that are well approximated by a 
modified polytrope of the form p = po + cP n . 

We find that highly detailed planet interior models, including temperature structure and phase changes, 
are not necessary to derive solid exoplanet bulk composition from mass and radius measurements. For 
solid exoplanets with no substantial atmosphere we have also found that: with 5% fractional uncertainty 
in planet mass and radius it is possible to distinguish among planets composed predominantly of iron or 
silicates or water ice but not more detailed compositions; with ~ 5% uncertainty water ice planets with 
> 25% water by mass may be identified; the minimum plausible planet size for a given mass is that of 
a pure iron planet; and carbon planet mass-radius relationships overlap with those of silicate and water 
planets due to similar zero-pressure densities and equations of state. We propose a definition of "super 
Earths" based on the clear distinction in radii between planets with significant gas envelopes and those 
without. 

Subject headings: extrasolar planets 



ner edge of its host star's habitable zone (jUdrv et alj 
I2007T ) . Microlensing surveys have discover ed two low-mass 
planets, ~ 5-5Mm ([Beaulieu et al l I2006D and ~ 13M e 
(|Gould et al.l l2006f ) at ~ 2.5 AU from their parent stars, 
suggesting that Neptune-mass planets are common. In one 
planetary system, the Neptu ne- mass planet is the most 
massive planet in the system (jGould et al.ll2006l h 

Space-based missions provide us with further motiva- 
tion for our study. COROlEI (ONES; launched 27 De- 
cember 2006) and Keplei0 (NASA; launch date 2008) will 
search for low-mass exoplanets that transit their host star. 
GALAQ (ESA; launch date 2011) will measure stellar dis- 
tances (and hence their radii) precisely, removing a limit- 
ing factor in deriving a precise planetary radius. Ground- 
based radial velocity techniques are pushing to higher pre- 
cision, and will enable mass measurements of many of the 
COROT and Kepler planets. 

To examine what can we determine about an exoplanet 's 
composition from measurements of its mass and radius, 
we derive theoretical mass-radius relationships for a wide 
range of exoplanet masses. To explore a wide range of 
masses and compositions, we make a major simplification: 
we make the approximation that the planet is at a uniform 
low temperature. This approximation of uniform low tem- 
perature serves as a practical simplification because the 
20051) and another with M p = 5M e orbiting at the in- equations of state (EOSs) are relatively well described at 
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1. INTRODUCTION 

The growing number and unexpected diversity of re- 
cently discovered extrasolar planets has motivated us to 
study the mass-radius relationship of solid exoplanets. 
The central question we pose is what can we determine 
about an exoplanet 's composition from its mass and ra- 
dius? The answer to this question requires numerical mod- 
els of planet interiors as well as an understanding of the 
current limitations and future prospects of the precision of 
planet mass and radius observations. 

The growing number of exoplanets includes some with 
interesting radii. The planet HD 149026 b has such a small 
radius for its measured mass that the planet must have a 
core with a mass of 60-70 M®, or 2/3 of the planet's to- 
tal mass (jSato et aT1l2005 | ). Another plane t, GJ 436b, is 
a Neptune-mass planet ([Butler et al.l 12004) that wa s re- 
cently discovered to show transits (jGillon et al.ll2007t ) and 
to have a Neptune-like radius (M p = 22.6 ± 1.9Mq and 
Rp = 3.95±0.35i?®). 

The unexpected diversity of exopla nets includes four- 
teen exoplanets with M v smi < 21Mm dButler et al. 20041: 
McArthur et al.l 12004 ISantos et all [20041: iRivera et al 



2005; 
2006; 



Bonfi ls et al l |2005[ iVogt et alj 120051: lLovis et al 



Udrv et all 120061 : iBonfils et all 120071: IMelo et al 



2007 ). including one with M p = 7.5M® ([Rivera et al 
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zero-temperature or at 300 K for a wide variety of materi- 
als below 200 GPa. The full temperature-dependent EOSs 
for the materials of interest are either unknown or highly 
uncertain at the temperatures massive solid planets can 
reach in their interiors and in the pressure range beyond 
the reach of static compression experiments (< 200 GPa) 
and the analytical high pressure laws of plasma physics 
(> 10 4 GPa). 

We adopt our approach f rom the foundational work 
of IZapolskv fc Salpeterl (|1969l ) who computed mass-radius 
relationships for homogeneous zero-t emperature spheres of 
single elements. We improve upon IZapolskv fc Salpeterl 
(1969) by using a more accurate E OS at pressures P < 
1000 G Pa. We further expand upon IZapolskv fc Salpeterl 
(|1969f ) by: considering more realistic planetary materials; 
considering differentiated planets; by exploring the effects 
of temperature on the planet mass and radius; and by in- 
vestigating potential observational uncertainties on planet 
mass and radius. 

This work complements the highly focused physical 
models of low- mass exoplanets offered by other authors. 
A detailed study by iLeger et alj (|2004f ) focused on water 
planets, and provided a detailed model of the interior and 
atmosphere of a 6 Mq planet with an interior composi- 
tion of: 3 Mq of water, 2 M® of a silicate mantle and 
1 M ffi iron and nickel core (see al so iKuchnerl (120031) for 
a description of a water planet.) IValencia et alj ((2006) 
calculated the mass-radius relationship for "super Earths" 
and "super Mercuries" (defined in their paper to be 1-10 
M s and similar composition to Earth and 1-10 Mercury 
masses and similar co mposition to Mercury, respectively) . 
IValencia et al.l (|2006l ) explore different mineral composi- 
tions of the mantle and core, and investigate whether the 
planets have solid or liquid c ores an d apply their model 
to Gl 876d in IValencia et al.l (l2007bl) and to deg e neracie s 



of planet interior comp osition in Valencia et al. (|2007af) 



Ehr enreich et alj (|2006l ) model small cold exoplanets to 
study the microlensi ng planet O GLE 2005-BLG-3 90Lb 
(|Beaulieu et al.ll2006l). In the pas t iSteyensonl d 19821) and 
more recentlv lFortnev et all (|2007l ) and lSotin et al.l (|2007t ) 
have also investigated mass-radius relationships of rocky 
and icy exoplanets. 

We take a broader view than these previous studies us- 
ing more approximate models in order to investigate plan- 
ets of a wide range of compositions and masses. We de- 
scribe our model in §2 and the equations of state (EOS) in 
§3. In §4 we present mass-radius relationships for homo- 
geneous and differentiated planets, and discuss the effects 
of phase changes and temperature. In §5 we describe a 
generic mass-radius relationship shared by all solid exo- 
planets under our approximation. We discuss the broad 
consequences of our study in §6 followed by a summary 
and conclusion in §7. 

2. MODEL 

We solve for m(r), the mass contained within radius, 
r, P{r), the pressure, and p(r), the density of a spherical 
planet from three equations: the mass of a spherical shell 

f =WKr); (1) 

the equation of hydrostatic equilibrium 
dP(r) ~Gm(r)p{r) 

dr r 2 



and the equation of state (EOS) 

P(r) = f(p(r),T(r)), 



(3) 



where / is a unique function for a given material. Dif- 
ferent approximations to the EOS have been derived and 
we describe our choice in detail in |j3j For the majority 
of the calculations for solid materials presented here, we 
neglect the temperature dependence of the EOS and use 
experimental data obtained at room temperature. The 
importance of thermal contributions to the pressure are 
discussed in tj4.2.21 

We numerically integrate equations (1) and (2) starting 
at the planet's center (r = 0) using the inner boundary 
condition M(0) = and P(0) = P CO ntrai, where P CO ntrai is 
a chosen central pressure. For the outer boundary condi- 
tion we use P(R P ) = 0. The choice of P C entral a t the inner 
boundary and the outer boundary condition P(R P ) = de- 
fines the planetary radius R p and total mass M p = m(R p ). 
Integrating equations (1) and (2) over and over for a range 
of ^central provides the mass-radius relationship for a given 
EOS. 

For differentiated planets containing more than one kind 
of material, we specify the desired fractional mass of the 
core and of each shell. We then integrate equations (1) and 
(2) as specified above, given a P CC ntrai and outer bound- 
ary condition. We switch from one material to the next 
where the desired fractional mass is reached, using a guess 
of the total planet mass. Since we do not know the total 
mass that a given integration will yield ahead of time, we 
generally need to iterate a few times in order to produce 
a model with the desired distribution of material. 

We tested ou r code by trying to duplica te the mass- 
radius icj : iryes_in Za^olsk^^SSi)itii] {l969), using their 
EOS (ISalpeter fc Zapolskvlll967ir Our mass-radius curves 
agreed with those in IZapolskv fc Salpeterl (|l969l ) to within 
a few percent. 

3. EQUATIONS OF STATE 

An equation of state (EOS) describes the relationship 
between density, pressure, and temperature for a given 
material in thermodynamic equilibrium. Because we com- 
pute models without temperature dependence, we choose 
a form of EOS that assumes uniform or zero temperature. 
For P < 200 GPa we use fits to experimental data, either 
the Vinet EOS or the Birch-Murnagham EOS (BME). For 
P > 10 4 GPa, where electron degeneracy pressure becomes 
increasingly important we use the Thomas-Fermi-Dirac 
theoretical EOS. In between these pressures the EOSs are 
not well known and we treat the EOS as described in : 



3.1. Low-Pressure EOSs: Vinet and Birch-Murnagham 

For pressures below approximately 200 GPa we rely on 
experimental data which have been fit to the com mon EOS 
formulae of either Vinet jVinet et aT]| 19871(19891 ) or BME 
(lBirchlll947t lPoirierl l2000l). For a derivation of these EOSs 
see lPoirierl (|2000D . The Vinet EOS is 



P = 3K or/ 



2/3 



exp 



r(^o-l) 



and the third order finite strain BME is 



(2) 



P = -Kn 

2 



4^0-4) 



- 1 



,(4) 



(5) 
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For the 4th order finite strain BME. 



-Kt 



»( 



the term 

2/3 _ f 



(6) 



x 



K K» + K>(K>-7) + ±f- 



is 



added to equation ©. Here r\ 
pression ratio with respect to the 
Kq = — V (§U) T is the bulk modulus of the material 



is the corn- 
density, 



= p/po 
ambient 

v \dv 

is the pressure derivative, and K(j is the second pressure 
derivative. The majority of experiments (from which Kq 
and K are derived) are typically limited to pressures less 
than 150 GPa and temperatures less than 2000 K. 

Both the BME and the Vinet EOSs are empirical fits 
to experimental data. The Vinet EOS is considered to be 
more suitable than the BME EOS for extrapolation to high 
pressures, because the BME is derived from an expansion 
of the elastic potential e nergy as a fu nction of pressure 
truncated to low orders (jPoirierl [2000). Where possible 
we choose the Vinet or BME fit provided for experimental 
data according to which fit best matches up with the TFD 
EOS at high pressures. In one case we used a fourth order 
BME where the term K' ' is determined theoretically (see 
gOD - Table 1 lists the K , K' , and the type of EOS fit 
we used for each material. 



3.2. High- Pressure EOS: Thomas-Fermi-Dirac 

The Thomas-Fermi-Dirac (TFD) theory was derived in 
the late 1920s as an approximate way to characterize the 
interactions of electrons and nuclei. The electrons are 
treated as a gas of noninteracting particles that obey the 
Pauli exclusion principle and move in the Coulomb field 
of the nuclei. Under the assumption that the potential is 
slowly varying, a self-consistent solution is derived so that 
Pauli exclusion pre ssure balances out the Coulomb forces 
(jEliezer et al.ll2002f) . 

These approximations lead to a comparatively simple 
description that works for any material and becomes in- 
creasingly accurate at high pressure. For each material, 
however, there is a pressure limit below which the TFD 
model is no longer valid, where the assumption of a nonin- 
teracting electron gas in a slowly varying potential breaks 
down. In real materials, the electrons occupy well-defined 
orbitals, which leads to chemical bonds and determines 
the crystal structure as a function of pressure and temper- 
ature. The TFD theory cannot describe chemical bonds 
and is insensitive to the arrangements of atoms in a par- 
ticular structure. At very high pressure, however, where 
the kinetic energy dominates over the Coulomb energy, all 
these effects become less important and TFD theory yields 
an increasingly accurate EOS. 

In this paper, we use a mod ified TFD theory developed 
bv lSalpeter fc Zapolskvi (jl967t ). The authors extended the 
original TFD theory by adding a density dependent cor- 
relation energy term that captures some of the interaction 
effects of the electrons. In all of the following TFD calcu- 
lations, we included the correlation energy correction cal- 
culate d with fit formulae provided in lSalpeter fc Zapolskvi 
(196TB We also follow this paper for the description of 
mixtures of different types of atoms. 



Since the TFD theory does not describe chemical bonds 
it does not reproduce the correct zero-pressure density, 
and can even be in er r or by up to a factor of two or more 
(|Zapolskv fc Salpeterl (|1969l ): see Figure [T]) . Instead we 
rely on experimental data for lower pressure, which are 
available for almost all materials we consider. 

3.3. Intermediate-Pressure EOS and Details for Specific 
Materials Used 

The pressure range from approximately 200 to 10 4 GPa 
is not easily accessible to experiment nor is it well de- 
scribed by the TFD EOS. Although shock experiments 
can reach pressures over 1000 GPa, a substantial contribu- 
tion to the pressure comes from the thermal pressure (see 
H4.2.2[) because the material is also heated under shock 
compression. When a shock wave passes through a solid 
material, its density is only increased by up to a factor of 
4. It is this limited compression ratio, which makes if very 
difficult to obtain low temperature EOS data beyond 200 
GPa, the range needed to model planetary interiors. For 
all materials but H2O in this pressure regime we simply 
use the Vinet/BME EOS up until the point where it in- 
tersects the TFD EOS curve, and the TFD EOS at higher 
pressures (see Figure [T]). A more accurate EOS in the in- 
termediate pressure range 200 to 10 4 GPa range demands 
new theoretical calculations. Such EOSs are not readily 
available since there are almost no applications that re- 
quire EOSs in this pressure range. 

Water Ice: As an example of how to fill this gap we 
used density functional theory to calculate the EOS of wa- 
ter ice in the phases VIII and X in the pressure range 2 
to 7700 GPa (Figured and Table 2). The theoretical EOS 
data presented her e are in agreement w ith water ice VII 
experimental data (jHemlev et al.| [l987) in the range 6 to 
127 GPa to within 3.5% in density for a given pressure. 
Our Vinet fit to the combined theoretical data for water 
ices VIII and X has parameters p — 1460 kg/m 3 , K = 
14.3771 GPa, and K' = 6.57972. This fit deviates from 
the tabular data in Table 2 by less than 2.5% in density 
for a given pressure. 

At P > 3 to 20 GPa water ice is in either phase VII or 
VIII, dependi ng on temperature (see t he wa ter phase dia- 
gram m, e.g., iPetrenko fc Whitworthl (jl999)). The struc- 
ture of water ice phases VII and VIII are extremely sim- 
ilar, differing only by the ordering of the hydrogen atom. 
This different structure causes a negligible difference in 
the EOS, making the experimental water ice VII and the 
theoretical VIII EOS comparable. At P = 60 GPa water 
ice VII/ VIII undergoes a phase change to water ice X. In 
this structure, the distinction between covalent bonds and 
hydrogen bonds goes away. Instead the hydrogen atom 
shifts to the mid-point between two oxygen atoms, while 
the hydrogen atoms occupy off-center sites in all ice struc- 
tures at lower pressures. More details on our choices of 
EOSs follow. 

For this work we adopt the f ollowing for our wate r EOS. 
We use the BME fit from the iHemlev et all (|l987l ) water 
ice VII data up to P = 44.3 GPa. At this pressure the the- 
oretical data and the experimental data agree precisely. 



■ as defined in S alpeter 



8 Note that the denominator in the second term of 

reproduces the EOSs in their paper. This definition of <j> differs from the one in 
factor of 1/3 1 / 3 . 



Zapolsky (1967) is listed as 4.3 1 / 3 ; the factor 4 X 3 1 / 3 actually 
Zapolskv fc Salpeterl l[l969), which appears to be missing a 
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Starting at P = 44.3 GPa we use the theoretical data 
derived from density functional theory, which represent 
state-of-the-art first-principles calculations. These calcu- 
lations were performed with the Vienna ab initio simula- 
tion package using the projec tor augmented-wave method 
(|Kresse fe Furthmullerlll996t ). The calculations predict a 
gradual transformation from ice VIII to X. The result- 
ing EOS is given in Table 2. At a pressure of 7686 GPa, 
our density functional theory calculation agrees with the 
TFD model and we use this pressure to switch to TFD for 
all higher pressures. In principle, density functional the- 
ory calculations can be performed for any material with 
a known crystal structure — this provides a way to bridge 
the gap in pressure between the experimental data and the 
TFD limit. 

For the liquid water EOS for P < 10 GP a, we use 
the logarithmic EOS (see, e.g., iPoirie n l200l. We use 



Kg = 2.28 GPa for seawater at 12° K from Hallidav et alj 
(120031) . 

Iron: For an Fe EOS we use the e phase of Fe with a 
Vinet fit up to P = 2.09 x 10 4 GPa. Note that the Vinet fit 
parameters are from ex perimental data with P < 330 GPa 
lAnderson et al. (2001). At this pressure the Vinet curve 
smoothly approaches the TFD EOS, and we switch to the 
TFD EOS. 

Silicate: For a silicate EOS we use the perovskite 
phase of MgSi03. We use a fourth order BME fit up to 
P = 1.35 x 10 4 GPa. At this pressure we switch to the 
TFD EOS. The fourth order BME fit is fro m a density 
functi o nal calculation up t o P = 150 GPa by iKarki et alj 
([20001) . IKarki et al.1 (|2000D note that their K' agrees with 
fits to experimental data from several sources and Kq is 
within the range of experimental data ( 247 GPa compared 
to 246-272 GPa). The advantage ofthe lKarki et alj (pOOO) 
fit parameters is that the fourth order BME is the only fit 
we found that smoothly matches the TFD EOS at high 
pressures. 

Other materials: Other m a terials used in this 
work include MgO (iDuffv et al.l fl995h. rMg.FelSiOs 
dKnittle fe Jeanlozl Il987h . and SiC (|Alexsandrov et alj 
1989). For these species we use the BME fit up to the 
pressure where they intersect the TFD curve. 

Carbon was the only material we used whose Vinet EOS 
and TFD EOS did not intersect; this was likely because 
we only considered the graphite phase at pressures below 
the TFD EOS. For graphite we interpolated between the 
Vinet and TFD EOS. For H, He and the carbon monoxide 
EOS, we used the TFD EOS at all pressures, for simplic- 
ity. We set the density to a constant at the low pressures 
(P < 10 8 Pa) where the TFD EOS is poorly behaved. 
We note that a C O EOS from shock experiments exists 
(jNellis et al.lll981l ) in the pressure range 5 to 60 GPa, and 
it has a density to within eight to sixteen percent of our 
TFD density. Our CO EOS and mass-radius relationship 
is therefore approximate. 

4. NUMERICAL RESULTS 

We now describe our numerical solutions to equations 
(1) through (3) using our assembled collection of EOSs. 
We used our model to investigate the mass-radius rela- 
tionships for planets from 0.01-1000 M ffi . The lower mass 
limit encompasses planets as small as Mercury and small 



bodies like the icy moons of Jupiter and Saturn. The 
upper-mass limit encompasses the 13 Jupiter-mass planet 
limit. Above this mass, self-gravitating H-He spheres un- 
dergo deuterium or sustained hydrogen fusion (depending 
on how massive the body is) and are not considered plan- 
ets. 

4.1. Mass-Radius Relationships 

4.1.1. Homogeneous Planets 

Building upon IZapolskv fe Sabeteri §969), we first con- 
sider planets of uniform composition. This artificial sce- 
nario helps us understand the fundamental properties of 
the planet mass-radius relationships. Figure 2] shows the 
mass-radius relationship for homogeneous planets of H, 
H/He (25% He by mass), H 2 (ice), MgSi0 3 (perovskite), 
and Fe. 

Homogeneous planets all show the same general trend in 
radius as a function of mass. For M p < 5OOM0 the plan- 
ets' radii increase with increasing mass. In this regime, 
Coulomb forces balance gravity in hydrostatic equilibrium. 
For large masses, M p 3> 5OOA/0, the compression in the 
interior is high enough to pressure ionize the atoms. At 
these large masses degeneracy pressure of free electrons 
balances gravity in hydrostatic equilibrium, an d as more 
mass is added to the planet, the planet shrinks ([Hubbard! 
Il984h . Although planets are not fully dege nerate (the term 
is res erved for stellar mass white dwarfs ([Chand rasckhar 
Il939f )). electron degeneracy pressure does have a signifi- 
cant effect on the mass-radius relationship for high plan- 
etary masses over 500 M®. In particular, planets of all 
compositions are approximately the same size for a decade 
of mass where the competing effects of Coulomb forces 

1/3 

(which cause R p ~ M p ) and electron degeneracy pres- 

— 1/3 

sure (Rr, ~ M v ) rough ly cancel each other out. See 
Zapolsk v fe Salpeterl (|1969f ) for a detailed discussion of the 
maximum radius for a given planet of homogeneous com- 
position. 

If we assume that our selection of materials spans all 
plausible major planet materials, then we can make some 
inferences from Figure [4] about the range of planet sizes. 
First, the Fe- planet mass-radius relationship shows the 
minimum radius a planet of a given mass can posses. Sec- 
ond, since water is the least dense of all the materials we 
studied (apart from H and He), the water planet curve 
in Figure [4] may serve to show the maximum radius for a 
planet with no substantial atmosphere. 

The mass-radius relationships for planets of homoge- 
neous compositions (Figure [3| can be used to infer the 
bulk composition of planets. Using the solar system as an 
example, and from Figures 2] and we could infer that 
Earth and Venus are composed primarily of a mixture of 
silicates and iron, while Mercury is composed predomi- 
nantly of iron. We could also infer that Uranus and Nep- 
tune are not giant H/He planets and nor are they "rock gi- 
ants" ; they are predominantly rocky or icy and must have 
small but significant gas envelopes. Jupiter and Saturn 
are grossly fit by the H/He curve, but the H/Hc interiors 
of hot Jupiters are dominated by thermal effects and are 
thus not fit well by cold homogeneous planets; indeed we 
are not aiming to model gas giant planets in this paper. 

Turning to exoplanets, Figure H] shows that the transit- 



5 



ing Saturn-mas^l exoplanet HD 149026b must contain a 
substantial fraction of elements heavier than H and He. 
More detailed evolutionary and interior models find that 
HD 149026b has 70 M m of roc ky material, almost 2/3 of 
its total mass (|Sato et alj|2005[ ). At 22.6 M® and 3.95 i?®, 
GJ 436b must have a significant H/Hc envelope because its 
radius is clearly larger than a pure water ice pl anet. The 
7.5 M® planet Gliese 876d (jRivera et al.ll2005h does not 
have a measured radius, but a radius measurement with 
a fractional radius uncertainty of 5% would distinguish 
among a predominantly rocky planet, a predominantly icy 
planet and a planet with a substantial gas envelope. 

We have modeled massive solid exoplanets up to 4000 
M® (up to 13 Jupiter masses). These planets would be 
Jupiter-mass planets composed of solid material. Such 
massive exoplanets are not yet known to exist. In the stan- 
dard planet formation theory massive planets are primar- 
ily composed of H/He and are limited to have rocky/icy 
cores of up to about 10M®. HD 149026b with a 70 M® 
core shows that a wider range of planets exist. Massive 
solid exoplanets of hundreds to thousands of Earth masses 
may be able to form around massive stars (B and O stars; 
5-120 M Q ) where the protoplanetary disk would contain 
enough heavy elements. Additionally these stars have high 
UV radiation and winds which could photoevaporate the 
nascent protoplanetary gas disk, allowing massive planets 
to form out of the remaining solid material. 

4.1.2. Differentiated Planets 

All solar-system planets have multiple layers of different 
compositions. These planets are differentiated, meaning 
the denser material lies beneath shells of progressively less 
dense material. We now consider differentiated planets of 
various compositions. We focus on materials that com- 
prise the solar system planets and moons: iron, silicates, 
H2O (ice), and H/He gas envelopes. We ignore elements 
that have abundances too low to affect our model for the 
planet radius. 

We explore two types of differentiated gas-free plan- 
ets, iron/silicate planets and water planets. These planets 
lack gas envelopes, although they may have atmospheres 
too small to affect the measured planet radius. Figure 0] 
shows the mass-radius relationship for differentiated plan- 
ets without gas envelopes. Figure [S] shows the same mass- 
radius relationships in more detail. The calculations as- 
sume a constant fractional mass in each layer. In general, 
the radii of differentiated planets (where the more dense 
components are interior to the less dense components) lie 
in between the radii of homogeneous spheres composed of 
the planet's most and least dense components. 

We investigate iron/silicate planets with iron cores and 
MgSiOa mantles. We consider Fe core mass fractions of 
32.5% ("super Earths") and 70% ("super Mercuries"). 
It is remarkable how well Mercury, Venus, and Earth's 
masses and radii are fit by these cold iron/silicate planets, 
as seen in Figure O We show the density as a function 
of radius for the silicate planets in Figure [5^ and b. As 
expected, the more massive planets have higher densities 
in their cores. 



We define water planets to be solid planets with > 25% 
H2O by mass. Jupiter's moons Ganymede and Callisto 
would be water planets by this definition^- We investi- 
gate water planets with iron cores and silicate mantles. 
We consider planets with fixed mass fractions: 45% wa- 
ter, 48.5% silicate s and 6.5% Fe (simil ar to Jupiter's icy 
moon Ganymede (|Schubert et al.l l2004j)): 75% water and 
22% silicates and 3% iron (the core and silicate shell mass 
ratio as Ganymede); and 25% water with a 58% silicate 
shell and a 17% iron core. Figure HJfc shows the density 
profiles as a function of planet radius for the 45% water 
planets. 

We compared our fiducial super-Earth model with a 
32.5% Fe core and 67.5% silicate mantle to a model that 
more closely represents Earth: a model with a 32.5% by 
mass core of FeS (where FeS includes 10% S by mass) 
and a mantle that includes 90% (Mg,Fe)Si0 3 and 10% 
MgO (|Poirierl [2000h . The results from this Earth model 
are shown by the squares in Figure [7| and agree closely 
with our fiducial super-Earth- like planet. It is remarkable 
how our simple fiducial model matches the Earth's radius 
to within three percent for a 1M S planet. 

There are further degeneracies among the mass-radius 
relationships for planets of different compositions. For ex- 
ample, planets with 10% water by mass (with 27% iron 
cores and 63% silicate mantles) have mass-radius curves 
that overlap with our silicate planet mass-radius curves. 
As a second e xample, if we ad opt a 6 M® water planet 
similar to the iLeger et al.l (|2004f ) water planet composed 
of 1 M® Fe, 2 M® silicates, and 3 M® H 2 0, we find the 
total planet radius differs by less than 0.5% from that 
of our model with 45% water, 48.5% silicates, and 6.5% 
Fe. While the total planet mass and radius are the same, 
the interior structure of the models are quite different, 
as shown in Figure [TJi (but note that the interior struc- 
tur e of our Leger-typ e planet is different than the one 
in (jLeger et al.l 12004) ). For the same mass, the radius 



of the Leger-type planet will be slightly lower, but both 
types of planet s fall alo ng the same mass-radius curve. See 
IValencia et al.l (|2007af ) and Li et al. (in prep) for detailed 
discussions on degeneracies in planets composed of iron 
cores, silicate mantles, and water outer layers. 

4.1.3. Planets with H/He Gas Envelopes 

We now turn to a discussion of differentiated planets 
with significant H/He envelopes. For simplicity and con- 
sistency, we use a zero-tem perature EOS for H and He 
(Sal peter fe Zapolskvi ri967). A zero-temperature EOS for 
a H/He mixture may represent real planets only poorly, 
but since a zero temperature EOS underestimates a gas's 
volume, using such an EOS allows us to one important 
point: adding a H/He shell can easily boost a planet's 
radius dramatically. For example, adding 20% H/He by 
mass can double a planet's radius. 

Figure [5] shows mass-radius relationships for planets 
with gas envelopes and fixed core masses. The cores con- 
tain 70% MgSiOa and 30% Fe. Five cases are shown: plan- 
ets with fixed core masses of 5, 10, 20, 50, and 100 M®. 
The values in Figure [8] are lower limits because thermal 



9 Saturn is 95 Mq and Jupiter is 318 Mgj. 

10 Under a water planet definition that includes H/He envelopes, Uranus and Neptune would also be water planets, because they are believed 



to have > 25% water by mass. 
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effects will move the gas curves up and left, i.e. to larger 
radii. For a detailed discussion of mass-ra dius curves for 
planet s with significant gas envelopes see iFortnev et alj 
( 2007]) and Adams et al. in preparation. In Figure [5] we 
also show the radius of the core of the H/He planets as 
a function of the total mass of the fixed-mass core and 
overlying H/Hc envelope. 

All planets with a significant amount of H/He will have 
a larger radius than the homogeneous water ice planets. 
Planets with even a relatively small fractional mass of 
H/He can therefore be distinguished from planets with an 
insignificant gas envelope (Figure [S| . Uranus and Nep- 
tune, for example, are believed to have about 10% of their 
mass in H/He material. Yet, Uranus and Neptune are 
up to two times larger than a planet of similar mass and 
core composition but without the H/He envelope. We pro- 
pose to call planets without a significant gas envelope — 
i.e. planets that lie below the pure water ice line — super 
Earths. 

4.1.4. Non-Standard Planets: Carbon and Helium 

Planets 

We now consider planets with compositions very dif- 
ferent from solar system planets, beginning with carbon 
planets. We have previously presented the idea of carbon 
planet s, planets composed of > 5 0% carbo n compounds by 
mass (jKuchner fc Seagerl ( |2006f ); See also ICameron et ail 
(|1988( ) and lGaidod (12000) ). Carbon planets should form 
in environments where the carbon-to-oxygen ratio C/O > 
1, in contrast to the solar abundance ratio (C/O = 0.5). 
When C/O > 1, the high-temperature condensates avail- 
able in chemical equilibri um are very different from the en- 
vironment with C/O < 1 (lLewislll974l:IWood fc Hashimoto! 
119931: lLodders fc Feglevl Il997f) . For example, SiC is the 
dominant form of Si instead of silicates (i.e., Si-0 com- 
pounds). Such carbon-rich environments may occur in a 
local area enriched in C or depleted in H2O in an oth- 
erwise solar-abundance protoplanetary disk. Carbon-rich 
environments would also occur in a protoplanetary disk 
with a global C/O > 1 like the disk that formed the plan - 
ets around pulsar PSR 1257+12 l|Wols zczan fe Frai]|[l 992'l. 
and do occur in the well-known Beta Pictoris debris disk 
(|Roberge et al.ll2006fl . 

The dominant composition of carbon planets is un- 
known. The details would depend on protoplanetary disk 
temperature, composition, relative abundance and depar- 
tures from chemical equilibrium. SiC and graphite are 
both plausible compositions (based on observations and 
calcul ations of carbon star atmospheres; lLodders fc Feglevl 
(|l997fi ). We compute the mass-radius relationship for 
three different kinds of carbon planets. The first carbon 
planet is a planet with an Fe core and an SiC mantle. The 
second is a carbon planet with an Fe core and a graphite 
mantle. The third type of carbon planet is a pure CO (car- 
bon monoxide) planet. A carbon monoxide planet could 
form in a stellar disk composed from a CO white dwarf 
that has be en shredded by its more massive stellar binary 
companion (|Livio et al.lll992ft . 

The main result from our carbon planet computations 
(Figure [5]) is that the carbon planet mass-radius relation- 
ships overlap those of the silicate and water planets. This 
is because the zero-pressure density of SiC (3.22 g/cm 3 ) is 



similar to that of MgSi03 (4.10 g/cm 3 ). While the zero- 
pressure density of graphite (2.25 g/cm 3 ) is almost twice 
that of water ice VII (1.46 g/cm 3 ), a graphite planet with 
an Fe core has a similar average density to a water planet 
with an iron core and silicate mantle. CO planets' mass- 
radius curves also overlap the mass-radius curves of water 
planets that contain iron and silicate because of a similar 
zero-pressure density. Our CO planets are similar in av- 
erage density to water ice planets, again because the zero 
pressure densities are similar. The precise mass-radius re- 
lationship of CO planets should involve temperature, and 
a more accurate EOS than our adopted one. 

We also compute the mass-radius relationship for a pure 
He planet. We again use t he zero-temperature He EOS 
(Salpctcr fc Zapolskvl I1967D to avoid introducing a free 
parameter based on the planet's unknown entropy. The 
helium mass-radius relationship is therefore approximate. 

A predominantly He planet may potentially form from a 
low-mass white dwarf. For example, a He planet can con- 
ceivably form in one type of symbiotic binary star called an 
AM CVn (AM Canes Venatici), composed of two very H- 
poor white dwarfs (i.e., He core white dwarfs) surrounded 
by a circumbinary helium accretion disk formed during 
mass transfer from t he less massive to the mor e massive 
white dwarf (see, e.g.. lPodsiadlowski et a l. (2003), and ref- 
erences therein). After it loses most of its mass, the less- 
massive WD may approach planetary mass. 

4.2. Phase Changes and Thermal Effects 

Real planets have phase changes and temperatures 
above 300 K in their interiors. Here we investigate the 
effect of phase changes and temperatures on the planet 
mass and radius. While models that include temperature 
and phase changes can be more realistic than our simple 
models, they complicate the interior boundary conditions 
and necessarily involve regimes wh e re the EOSs ar e highl y 
uncertain. Se elLeger et all (120041) ; IValencia et"al] (120 06); 



iSelsis et all (|2007t ): ISotin et all (|2007t k IValencia et al l 
(2007b?) for low- mass planet interior models that do in- 
clude phase and temperature changes for specific types of 
planets. 

4.2.1. Phase Changes 

In this section we show that the low-pressure phase 
changes have little effect on the planet mass-radius curves. 
This lack of effect is because the low-pressure phase 
changes occur at pressures < 10 GPa. Figure [TT1 shows the 
pressure at m{r) = 0.97Af p for homogeneous planets. Fig- 
ure QT] shows that for planets > 3Mq, most of the planet's 
mass is at pressures higher than 3 GPa (with the exception 
of pure water planets) . Here we discuss the phase changes 
in some of the key planetary materials individually. 

Water: We begin by investigating the liquid water to 
water ice phase change. The room-temperature bulk mod- 
ulus Kq for liquid water is an order of magnitude smaller 
than Kq for water ice VII. Furthermore, the densities of 
liquid water and water ice VII are different by 50% (see 
Table 1). Despite these large differences in physical prop- 
erties, it is important to remember that liquids are still 
not highly compressible materials and we expect them to 
behave more similarly to solids than to gases under high 
pressure. In order to investigate the effect of a liquid wa- 
ter phase for a water planet, we consider liquid water to 
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be present at P < 10 GPa. The temperature would be 
~ 650 K for water to remain a liquid at this high pres- 
sure. We consider the differences between a pure water 
ice planet with and without the liquid ocean. We find less 
than three percent difference in the radius for water ice 
planets without and with a liquid water ocean for planets 
in the mass range 1 to 4 Ai®. We find < 1% difference in 
the radii of the two types of planets for M p above 5 M® 
(see Figure [7]). 

Silicates: We now turn to phase changes in our MgSiC>3 
perovskite planets by considering the low-pressure phase 
of MgSiC>3 called enstatite. We adopt the low-pressure 
enstatite phase at pressures less than 10 GPa. The differ- 
ence in Kq and po between the high-pressure perovskite 
phase and the low-pressure enstatite phase are much less 
than the differences between liquid water and water ice 
(see Table 1). Hence, just as in the water planet case, 
we see a relatively small difference (< 1%) in radius be- 
tween silicate perovskite planets with and without a layer 
of enstatite at P < 10 GPa (see Figure CJ). 

Iron: We now turn to phase changes in Fe. The low- 
pressure phase of Fe is Fe(a). Fe(e) and Fe(a) have the 
same Kq and K to within the experimental error bars, and 
only slightly different zero-pressure densities (see Table 1). 
Additionally, Fe (a) only exists at < 10 GPa; whereas in 
a differentiated planet Fe is expected to exist at pressures 
higher than 10 GPa and in a homogeneous Fe planet most 
of the planet's mass is above 10 GPa (Figures [TtJl and |TT|) . 
For these reasons the phase change of Fe(a) to Fe(e) has 
little to no effect on the mass-radius relationship of Fe 
planets. 

High pressure phase changes: Even if there is a 
phase change at high pressures (not considered here), we 
expect the associated correction to the EOS to be small, 
and hence the derived planet radii to be reasonably accu- 
rate. Phase changes arise from a rearrangement of atoms 
in the crystal structure and the associated modifications 
in chemical bonds. With increasing pressure, the atoms 
become packed more and more efficiently and the im- 
portance of chemical bonding patterns drops significantly. 
For example, above 10 GPa, the postperovskite phase 
(|Murakami et al.1 |2004[ ) and perovskite phase of silicate 
have very similar zero pressure vol umes and bulk modul i, 
to less than a few percent difference (|Tsuchiva et a l. 2005) . 
We caution that although phase changes do not need to be 
considered, theoretical calculations are needed to compute 
an accurate EOS at high pressures (i.e., between 200 and 
10 4 GPa) because in many cases an extrapolation of the 
BME or Vinet EOS fit will not do (see g3T3j) . For exam- 
ple, for some water ice planets with masses above 10 M®, 
the effect of extrapolating the BME fit without consider- 
ing a more accurate higher pressure EOS is considerable 
(Figure EH). 

At very high pressures the TFD EOS becomes valid. 
To understand the planet mass range where TFD effects 
are important, we show planet mass-radius relationships 
in the case of ignoring the TFD EOS in the dotted lines 
in Figure [JJ i.e. in using only the extrapolated BME or 
Vinet fit. 



4.2.2. Temperature 



Temperature has little effect on the radius of solid exo- 
planets. The reason is that the density of a solid changes 
by a relatively small amount under the influence of thermal 
pressure. At low pressures, the crystal lattice structure 
dominates the material's density and the thermal vibra- 
tion contribution to the density is small in comparison. 
At high pressures, the close-packed nature of the atoms 
prevents any significant structural changes from thermal 
pressure contributions. Moreover, any change in average 
density of a planet of fixed composition results in a smaller 
change in the planet radius than the change in average 
density, because R p ~ pp -1 / 3 . 

Some authors have stated or shown with models that 
temperature is not important for deriving a p l anet's 



total mass and radius (see 



Valencia et at! 120061 



I2007bt iFortnev et all l2007t ISotin et alj 120071 and refer- 
ences therein). In this subsection we go into detail to 
estimate the effect of temperature on density for the three 
main materials that we studied, iron, silicates, and water 
ice VII. 

Heating planetary material adds a thermal contribution 
to the pressure term. In a first approximation one can as- 
sume that the thermal pressure is a linear function o f tem- 
perature fe.g.. lAnderson fe GotolH989tlPoiriejl2000l^ . The 
thermal contribution to the pressure, Pth, can be written 
as: 

Pth dT = £ a(T)K (T)dT. (7) 

Here a(T) is the thermal expansivity and Kq{T) is the 
bulk modulus as before. 

For most cases the term a(T)K Q (T) is indepen- 
dent of volume above the Debye temperature 0_d 
(|Anderson fe Gotolll989f ). The thermal pressure term can 
then be written 

r&o 

P th = aK dT + aK (T-e D ). (8) 



Note that the 0_d is around 700K for silicate- and Fe- 
bearing minerals. Od is significantly lower (300K) for wa- 
ter ice and significantly higher for carbon (2000K). More- 
over, for Earth-bearing minerals I Anderson fe Go to (1989) 
find that Pth is linear in T down to much lower tempera- 
tures, and find that 

P th = aK Q (T~ 300 K). (9) 
For (Mg,Fe)SiQ 3 we use a K„ = 0.00692 GPa/K from 
lAnderson fe Masudal (|1994f ). We note that, for our pur- 
poses, (Mg,Fe)Si03 is similar enough to MgSiOs for esti- 
mating the thermal pressure. 

For metals the thermal excitation of electrons must 
be taken into account b y including a higher order term 
(Isa ak fe And erson 200^, and references therein), 

P th = aK (T ~ 300 K) + (T - 300 K) 2 . (10) 



We use values of aKn 



dT v y 
= 0.00121 



GPa K" 



and (daKo)/(dT) y = 7.8 x 10~ 7 GPa K~ 2 from 



anu \uai\Qi/ \ui }y = 

llsaak fe Andersonl (|2003l ). 



For H 2 ice VII, we constructed P t h based on (P,V) 
isotherms according to t he therm o dynam ic relations in 
equations (2) and (3) in iFei et al.l (| 19931 ). This method 
uses a linear fit with T to a(T) and a parameter ryi 
which describes the pressure effect on the measured vol- 
ume at high temp erature. We take the parameters from 
iFrank et al.l (|2004l ). 
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Figure [T2] shows Pth vs. temperature for Fe, H2O ice 
VII, and (Mg,Fe)SiC>3. P t h increases more slowly at high 
pressure than at low pressure because at high pressure the 
thermal expansion becomes a constant; the atoms become 
tightly packed so that any thermal pressure has a decreas- 
ing contribution to the total pressure. Although water ice 
da ta is only availabl e for ice VII to 50 GPa and ~ 800 
K (|Frank et al.ll2004T ). we expect it to show the same Pth 
trend with increasing temperature as other materials. 

We can estimate the change in density caused by a ther- 
mal pressure. We computed p(T) corresponding to the to- 
tal pressure Ptotai = P + Pth using the EOSs described in 
13.31 Figures [T2] to [T3] show (P,p) isotherms for the main 
materials we studied. 

For MgSi03, the decrease in density due to thermal pres- 
sure is less than 4 percent over the temperature range 
300 K to 6000 K and above 10 GPa where the material 
is solid (Figure flUa). Below 10 GPa the decrease in den- 
sity is less than 2.5 percent for temperatures up to 1200 K 
(Figure [T3b and for densiti es less than the zero-pressu re 
300 K density see Figure 7 in lAnderson fe Masudal ((1994).) 

Fe is known to have a higher thermal expansivity than 
MgSiO-3 and its density therefore decreases more than 
MgSiOa for the same temperature increase. The Fe den- 
sity decrease is less than four percent for pressures above 
100 GPa (Figure fT4k). Earth's Fe core is at pres sures 
higher than 100 GPa (jDziewonski fe Andersonlll981f ). and 
we expect the same for more massi ve differentiated plane ts 
with similar composition (see, e.g. IValencia et al.l (|2006[ )1. 
We further note that as the planet's Fe mass fraction in- 
creases, more of the planet's mass is at higher pressures, 
making the material's density change at low pressure in 
response to temperature less relevant. For example, Fig- 
ure [TT] shows that for a pure Fe planet with M p > O.IM®, 
97 percent of the planet's mass is above 1 GPa and for 
planets with M p > 1M® 97 percent of the planet's mass 
is above 10 GPa. 

Us ing available data to 50 GPa and 800 K ([Frank et alj 
2004), we find that H2O water ice VII density changes 
by less than a few percent throughout this pressure and 
temperature range (Figure [Hp). We note that above ap- 
proximately 1000 K wa ter ice reaches an ionic fluid phase 
(|Goncharov et al1l2005T ). Nevertheless, we emphasize that 
our water ice EOS (which includes phases VII, VIII, and 
X; see Table 2) agrees with recent Hugoniot shock data to 
within the experimental uncertainties. This recent data 
is from iLee et ail (|2006f) and ranges in pressure from 47 
to 250 GPa and temperature from 2100 K all the way to 
19,000 K. More work needs to be done to quantify the 
thermal pressure effects above 250 GPa and in the ionic 
phase, which is beyond the scope of this paper. 

Our (P, p) isotherms in Figures Q2] and [H] can be in- 
terpreted in light of the temperature-pressure profiles of 
planetary interiors. Earth, for example, has temperatures 
of approximately 1000 K at 5 GPa 2000 K at 150 GPa 
and up to 6000 K at 350 GP a (lPoirierll2000D Supe r-Earth 
interior models calculated in IValencia et al.l (|2006T ) shows 
that below 1 GPa the temperature is less than 1600 K, 
above 100 GPa the temperature is 2500 K, and the tem- 
perature is up to about 7000 K at hundreds of GPa. Isen- 
tropes for Neptune are expected to be about 2500 K at 20 
GPa and to reach 7000 K at 600 GPa. The fact that high 



temperatures are reached only at high pressures lessens 
the temperature effect because the fractional contribution 
of the thermal pressure to the total pressure decreases as 
the total pressure increases (Figure [12). 

We argue that even for the short-period hot 
exoplanets — those in few day orbits whose surface tem- 
peratures could reach 1000 to 2000 K — the thermal pres- 
sure contribution to radius is still small. This statement is 
based on the argument that for solid planets above one to 
a few Earth masses (depending on composition) the frac- 
tional radius affected by such high temperatures is small. 

Based on the above discussion of density decrease as a 
function of temperature we illustrate the effect of temper- 
ature on planet radius by the following example. Taking a 
case where the average density is overestimated by 3.5 per- 
cent in a uniform 300 K temperature models, the planet's 
total radius would be underestimated by only 1.2 percent. 
This is due to the R p ~ Pp^ 1 ^ 3 scaling. 

5. A GENERIC MASS-RADIUS RELATIONSHIP FOR SOLID 
EXOPLANETS 

A glance at Figure 2] may suggest to the astute reader 
that the mass-radius relationships for a variety of plan- 
ets all have a similar functional form, perhaps caused by 
some symmetry of the underlying equations. Indeed, there 
is such a symmetry, and a common functional form: a 
generic mass-radius relationship which we describe here. 
This generic mass-radius relationship is valid for planet 
masses up to about 20 M®. 

5.1. A Modified Polytropic Equation of State 

Our generic mass-radius relationship is based on the 
similar forms of the EOSs of all solid materials we have 
considered. The zero-temperature or 300 K temperature 
EOSs for the solid materials we considered can be approx- 
imated by the function 

p(P)=p + cP n . (11) 

Table 4 lists the best-fit parameters po, c, and n for some 
materials over the range P < 10 16 Pa. The density given 
by these approximate EOSs match the density given by the 
more detailed EOSs we used above to within 2 to 5% for 
the pressure ranges P < 5 x 10 9 and P > 3 x 10 13 Pa. At 
intermediate pressures, the discrepancy ranged from less 
than 1% to 12%. 

The similarity of the EOSs of all solid materials stems 
from the the behavior of chemical bonds under pressure. 
At low pressures the chemical bonds of the material can 
withstand compression. Above some pressure, the energy 
imparted to the material breaks the chemical bonds and 
the material structure radically changes. The "crossover" 
pressure is roughly the material's bulk modulus. For our 
simple analytic EOS, the bulk modulus is 



K 



(12) 



Our approximate EOS is a modified polytropic EOS. 
Polytropic EOSs are of the form P = K p p 1+1 ^ m , or 
p = (P/K p ) m /( m+1 \ where K p is a constant and m is the 
polytropic index. Our approximate EOS fit differs from 
the polytropic EOS by the addition of the constant po; 
it incorporates the approximate incompressibility of solids 
and liquids at low pressures. 
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5.2. Dimensionless Equations of Planetary Structure 

The new generic EOS contains some dimensional quan- 
tities, p and c, that allow us to conveniently write equa- 
tions (1) and (2) in dimensionless form. We rescale the 
variables P, p, m, and r as follows, where the subscript s 
refers to a scaled variable, 

P. = - (13) 
Po 

Ps = ^r (14) 



Pi 

r 

rn 
mi 



Pi 



1/n 



(15) 
(16) 

(17) 
(18) 

(19) 

With this change of variables, the equations of planetary 
structure become: mass of a spherical shell 

dm s (r) 



where 



and 



n = G-^p^- 1 ^- 1 ^ 



mi = r 1 p . 



dr s 

hydrostatic equilibrium 

dP s {r s ) p s (r s )m s (r s ) 



dr. 



and the EOS 



^rip s {r a ); (20) 

(21) 

p s = l + P s n . (22) 

Figure [15] shows dimensionless mass-radius relationships 
derived by numerically solving the above equations the 
same way we solved the unsealed equations. It shows 
the total dimensionless planet mass, M s as a function of 
the scaled planet radius, i? s , where R s is defined by the 
outer boundary condition P(R S ) = 0. At R s we also have 
M s = m s (R s ). We generated the numerical mass-radius 
curves by solving the equations for a range of central pres- 
sure values of P s . 

The dimensionless mass-radius curves depend only on 
n. We plot curves for three values of n: n = 0.513 (H2O), 
n = 0.528 (Fe), and n = 0.544 (silicate). These values of 
n span the range of behaviors of all the EOSs we stud- 
ied. We also show solutions for the interior structure of 
homogeneous planets in Figure [TBI 

The solutions behave quite differently on either side of 
the line M s = 1. For M s < 1, R s strictly increases with 
M s , and does not depend on n. For M s > 1, R s depends 
strongly on n and does not necessarily strictly increase 
with M s . This contrasting behavior arises because for 
M s <C 1, the EOS reduces to p s = 1, and for for M s > 1, 
the EOS reduces to the polytropic form p s — P™. 

Table 4 lists some values of mi, ri, and Pi for some 
EOSs based on the values of po, n, and c listed in Ta- 
ble 3. Curiously, the scaling parameters mi and n are 
somewhat similar for the polytropic-like mass-radius solu- 
tions for H2O and Fe because these materials have similar 
ratios of K /p . 



For M s < 4, the dimensionless mass-radius relationship 
is approximately 

logi R s = ki + l/31og 10 (M s ) - k 2 M k s \ (23) 

where h = -0.20945, k 2 = 0.0804, and k 3 = 0.394; 
this approximation is good to 1% over this range. For 
M s > 4, the scaled radius becomes a strong function of 
n. But We can use the analytic function in equation ([2"3"| 
to describe the dimensionless mass-radius relationships for 
scaled masses, M s , up to ~ 40 if we use the values of fcj 
for different materials listed in Table 4. 

Equation ([23]) and the scaling parameters in Table 4 
provide a convenient approximate summary of the results 
of this paper for homogeneous planets. 

It may appear that the scaled mass-radius relationship 
is not useful for differentiated planets, since differentiated 
planets combine materials with different equations of state. 
We find, however, that even the radii of differentiated plan- 
ets are well described by the functional form shown in Fig- 
ure [T5] for M s < 4 (i.e., up to where the scaled numerical 
solutions for different materials differ from each other). 

This circumstance provides us with a convenient way to 
summarize our results for differentiated planets. To any 
differentiated planet model, we can assign an effective mi 
and rx. These effective scalings allow us to summarize 
all of our calculated mass-radius relationships using equa- 
tion ([2"3"]l — regardless of whether they were computed with 
the modified-polytropic EOS or our actual Vinet/BME 
and TFD EOSs. Table 4 lists some effective values of mi 
and r% for a few examples of differentiated planets. We 
calculated these numbers by comparing the mass-radius 
curves for differentiated planets (shown in Figure [3} to 
equation (|23[) . To find the approximate mass-radius rela- 
tionship for any given planet, look up mi and r% in Ta- 
ble 4 and plug these numbers into equation ([23]) using 
h = -0.20945, k 2 = 0.0804, and k 3 = 0.394. 

The scaled variables help us to understand why the solid 
planet mass-radius relationships are very nearly the same 
for all materials we considered, for M s <= 1. We first re- 
call that M s = 1 is defined as the mass so that the central, 
i.e. the maximum, pressure in the planet is P s = 1. Next, 
for P s < = 1, i.e. everywhere in the planet for M s <= 1, 
the density (or the scaled density) never changes by more 
than 2.5 percent if a variation of n in the range of 0.513 
to 0.549. Hence the radius (or the scaled radius) never 
changes by more than 0.85 percent (due to the R p ~ p -1 / 3 
scaling) . 

5.3. Analytic Treatment of the Dimensionless Equations 

Here we will derive an approximate analytic solution to 
the dimensionless equations of planetary structure. The 
existence and form of the solution demonstrates why the 
mass-radius curves for various planets all look so similar. 
The good agreement between this approximate solution 
and our calculations gives us confidence in our results. 

We first discuss a general analytic solution to equations 
(1) and (2), followed by an application to the dimension- 
less equations (|2T))) and ([2"Tj) . For most equations of state, 
equations (1) and (2) cannot be solved analytically, even 
given the approximation of zero temperature. We can in- 
corporate two ideas to allow new analytic progress. First, 
over a wide range of low pressures below a GPa, solids and 
liquids change their densities by a small amount much 
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less than ten percent. This point enables us to assume an 
equation of state that incapsulates the idea of materials 
that are largely incompressible over a wide range of low 
pressures: 

p(P)=p + f(P). (24) 

We have shown in H5. II that f(P) = cP n is a good approx- 
imation. We further note, however, that in the range of 
low pressures we can assume f(P) <C po- 

The second point that enables an analytic treatment is 
that, when a planet is massive enough that it begins to 
compress under its own gravity, the compression is most 
acute at the planet's center. With these two ideas in mind, 
we can obtain an approximate solution to equations ([I]) 
and (J2), and a mass-radius relationship for the case of a 
low-mass, slightly compressed planet. 

To zeroth order, p(P) does not depend on P, so we can 
integrate equations fl} and @ to find the zeroth order 
solutions: 

m(r) = -nr 3 p , (25) 



Here, 



P(r) « Pc - ^Gr 2 p 2 . 



P c --nGR pPo -—^, 



(26) 
(27) 



where M p is the total mass of the planet, and R p is the 
planet's radius. 

Now we can use the zeroth order solution for P(r) to 
write equation l|24p to 1st order in the term /(P(r)): 



p(P) *p Q + f(P(r) = P c 



-nGr 2 p 2 ). 



(28) 



Keeping in mind that the compression is most important 
in the center, we will expand this expression for p about 
r = 0, to get an expression for the density accurate to first 
order in /(P(r)) and second order in r/R: 



P(Pc 



Po 



f(Pc 



nGr 2 plf\P c ) 



(29) 



If the compressed region is confined to small r/R, then we 
can use this expression as an approximation for the pres- 
sure everywhere in the planet. We can then substitute this 
expression into equation |l} and integrate to get 



m(r) 



p(P c ) - -7rGr 2 p 2 f(P c ) 



(30) 



When we evaluate this equation at r 
desired mass-radius relationship 



R p , we find the 



p(Pc) - \nGR 2 pP lf'(P c ) 



The mean density is 

P = P(Pc] 



-nGR 2 p p 2 Q f'(P c ). 



(31) 



(32) 



We can see that p(P c ) > p > Po, i-e., the mean density 
is higher than the zero pressure density but lower than 
the density near r — 0. Also, for high P c the radius will 
decrease (this is the case for f(P c ) = cP™, the modified- 
polytropic form in equation (JTTJ)) . 



We can evaluate equation (|52")l to first order in /(P(r)) 
using equation (f27|) . In other words, the mean density of 
the planet is approximately 



p = p(P c )--f(P c )P c , 



(33) 



where P c is given in equation (|27p . 

We now apply the above analytic equations equa- 
tions (f2"6")) - (P3")) to our dimensionless equations. We begin 
by considering our scaled equation of state as 



Ps (p s ) = i + p: 



(34) 



We proceed under the assumption that P s <C 1 and there- 
fore P™ <C 1. For the dimensionless equations we find: the 
scaled pressure 



where 



Ps(r s ) «P, 



P.,c = gTTift 



the scaled central density 

p S (P S ,c) « 1 



the scaled mass 



m s (r s ) ps -Trr 3 s 



1 + P™ 



2 

— imr 

5 



2pn- 



the desired mass-radius relationship 



M, 



-TiRi 



1 



and the average density 

Ps = 1 - 



(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



Because we chose P s <C 1 for this analytic derivation, 
we know exactly over what range of parameters these ap- 
proximations are valid. They apply where P 8 < 1, which 
is also where R s <C 1 and M s <C 1. The correction term 
in equations (J39J) and gO]), [(l - |n) (1^)™] < 1. The 
fi s <C 1 limit therefore shows why the scaled mass-radius 
relation depends very weakly on composition: the correc- 
tion term in equation (|39[) is small. 

If we consider the scaled mass-radius relationship equa- 
tion (f3"9")l slightly beyond where it is formally valid (R s < 1 
instead of R s 1), we find that it is still a reasonable ap- 
proximation. While the scaled mass-radius relationship 
equation (|3"9")) is good to within 1% at M s = 0.245 it is 
good to within 5% at M s = 0.36 (compared to the numer- 
ical solution to the scaled equations). Even at R s ps 1, 
the correction term is never larger than (1 — |n)(|7r) n . 
Since n ranges typically from 0.5-0.6, this maximum value 
for the correction term ranges from about 0.997 to 1.013. 
In other words, even for R s ps 1, the correction term in 
equation ([39]) varies over only a range of about 0.02 as a 
function of planet chemistry. 

Figure [T5] compares this analytic mass-radius approxi- 
mation to the full numerical solution of the scaled equa- 
tions, justifying that our calculations are correct. 
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6. DISCUSSION 

6.1. Exoplanet Mass and Radius Observational 
Uncertainties 

We now discuss the observational uncertainties on the 
mass and radius of transiting exoplanets. We adopt a con- 
servative range for planet mass and radius fractional un- 
certainty of 2-10%. The 10% limit is the uncertainty range 
for typical exoplanet detections, arising from measurement 
uncertainty due to the data quality. Once discovered, most 
exoplanets of interest will be followed up with larger tele- 
scopes and/or more observations to refine the planet mass 
and radius beyond that determined from the observational 
discovery. If measurement uncertainty is not a dominant 
factor for planet mass and radius uncertainty, then the 
stellar mass and radius uncertainty are. This is because 
the exoplanet mass and radius are derived from quantities 
that involve planet-star mass or radius ratios. Note that 
assuming the stellar noise to be small, the measurement 
uncertainty and stellar mass or radius uncertainty add in 
quadrature. 

Our estimate of the 2% planet mass and radius uncer- 
tainty is based on an optimistic assessment of the 2% stel- 
lar mass and radius uncertainty likely to be possible for 
millions of stars in the future. This kind of high-precision 
measurement will be enabled with accurate distances and 
precise stellar fluxes by the GAIA space mission (ESA; 
launch date 2013). In practice the radius can be inferred 
directly from the stellar fluxes and distances; in principle 
a precise radius is limited by the correction from a mea- 
sured stellar flux to the star's bolometric flux (D. Sasselov, 
2006 private communication). In contrast to our 2% best 
case scenario, current typical stellar ma ss and radius un- 
certainties are on the order of 5-10% (iFqrd et ail [l999; 
ICodv fc Sasseiovl 120021 : iFischer k Valentil I2005T) . These 
stellar masses and radii are derived from in terior and evo- 
lution model fits to observed stellar spectra. ISozzetti et al.l 
(|2007ft show that a more precise stellar radius (and hence 
planet radius) can be derived using stellar evolution mod- 
els in in the a/i?» vs. Tejf parameter space instead of 
in the usual log g vs. T e ff parameter space. This is 
partly because a/R* (a measure of stellar density) can 
be determined with high precision from the planet tran- 
sit light curve for planet's with zero orbital eccentricity 
(Sca ger k Mallen-Orn elas 2003}). Provided the photome- 
try is good enough, few percent uncertainty in star and 
planet radii may become routine. We note that a dif- 
ferent technique, interferometry, can measure stellar radii 
directly, but the current uncertainties are much higher 
than 2% and the technique is limited to a small number 
of nearby stars. 

The mass determination of low-mass planets by ground- 
based radial velocity techniques require hundreds of ob- 
servations. Mass deter mination for many exoplanets may 
therefore be inhibited (|Lovis et alj 1200a ) until a number 
of dedicated ground-based telescopes are available. Given 
such limitations of telescope time and current technology, 
the optimistic future exoplanet mass uncertainty range is 
likely closer to 5-10% for the low-mass solid planets of 
interest (< 20M ffi ). With current technology and due to 
the faintness of the host stars, many low-mass exoplan- 
ets discovered from transit surveys (e.g., COROT or Ke- 



pler) will not have measured masses at all. As an example 
of what it would take to detect an Earth-mass planet in 
a 50-day orbit about one of the brightest sun-like stars: 
one eight meter diameter telescope dedicated to five bright 
stars monitored over five years (R. P. Butler 2006, private 
communication). Earth-mass planets more distant from 
their host stars, such as Earth-like orbits about sun-like 
stars cannot be detected with any current technology. 

6.2. Possible Exoplanet Compositional Distinctions 

Given the observational uncertainties, what composi- 
tional distinctions among exoplanets can we make? With 
an upper limit of 20% on the planet mass uncertainty we 
will be able to say robustly that the planet is predomi- 
nantly composed of solids or if it instead has a significant 
gas envelope as do Uranus and Neptune. 

The following discussion assumes we can ascertain that 
the planet has no substantial atmosphere or envelope that 
would contribute to the planet radius (c.f. Adams et al., 
in prep). 

With a 10% uncertainty we may be able to comment 
on the presence of a large amount of water or iron — if the 
planet fortuitously has a very low density or a very high 
density within the radius range for solid planets. 

A ~ 5% uncertainty will allow us to distinguish 
among planets composed predominantly of water-ice, pre- 
dominantly of silicates, and predominantly of iron (sec- 
tions 14.1.11 and I4.1.2p . These planets are relatively well 
separated on the mass-radius diagram (Figure [4]) . Their 
separation is primarily due to the low density of water ice, 
the intermediate density of silicates, and the high density 
of Fe. Even with 5% exoplanet mass and radius uncertain- 
ties, it is not possible to identify the detailed composition 
such as the fraction of different material in the core and 
differentiated layers. 

Identification of water planets is possible with ~ 5% 
planet mass and radius uncertainty if the planet has more 
than 25% water ice by mass and plausible iron-to-silicate 
ratios. With the same ~ 5% uncertainty, water planets 
with 50% water by mass with any iron-to-silicate ratio can 
be identified. Water exoplanets should exist; in our own 
solar system Jupiter's satellite Ganymede is 45% water ice 
by mass. Water ice planets should be the easiest of the 
solid exoplanets to detect observationally because of their 
large radius for a given planet mass — they c ould be as large 
as 3Rm for M p = 2OM . We note that IValencia et alj 
( 20075 came to the same conclusion for the one planet 
mass they explored water compositions of, the 7.5 Mq Gl 
876d. A detection of a low-density water-ice planet orbit- 
ing far interior to the expected snow line in the disk where 
the planet formed would be strong evidence for planet mi- 
gration. 

With 2% uncertainty in planet mass and radius we 
would be able to determine not only the basic composi- 
tion of any iron/silicate/water planet, but constrain the 
relative fraction of each material. However, even with 
2% uncertainty in planet mass and radius many degen- 
eracies in planet composition remain for a planet of the 
same mass and radius. For example, carbon planets with 
silicate mantles and iron cores are indistinguishable from 
silicate planets with small iron cores, or with planets that 
have a small (~ 10% by mass) fraction of water on top 
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of a silicate mantle and iron core. As a second example, 
planets with deep water oceans (even 100 km deep) on any 
kind of planet are not identifiable with even 2% fractional 
uncertainty in the radius. This is because a liquid water 
ocean contributes only a small amount to the overall planet 
radius compared to a pure water ice layer. For a third ex- 
ample, see Figure BUI and the accompanying discussion. 
See also I Valencia et all (|2007a[) for a detailed discussion of 
degeneracies in composition for iron-silicate-water planets 
of 1-10 M m . 

Observational uncertainties in planet mass and radius 
are unlikely to be better than a few percent in the next 
decade. We therefore argue that detailed exoplanet inte- 
rior models are not needed to infer exoplanet bulk compo- 
sition. The analytical form derived in this paper (equa- 
tion with scaling relations provided from model 
curves (e.g., in Table 4), should be sufficient for the near 
future. 

Exoplanet atmosphere measurements may inform us 
about the planet interior, removing some of the degenera- 
cies in a given solid exoplanet mass and radius relation- 
ship. Therefore, despite our conclusions that detailed inte- 
rior planet models are not needed to determine exoplanet 
bulk composition, we point out that detailed interior mod- 
els are required to understand planetary atmospheres. In 
turn, observations of exoplanet atmospheres may help us 
infer more about the exoplanet interior composition. At- 
mosphere measurements will also help us to identify in- 
teresting features on planets, such as the presence of a 
deep liquid ocean. A saturated water vapor atmosphere 
would be a better indicator than an exoplanet mass and 
radius, considering we cannot identify deep water oceans 
with mass-radius relationships. 

7. SUMMARY AND CONCLUSIONS 

We have modeled cold solid planets of a variety of com- 
positions including iron, silicates, water, and carbon com- 
pounds. The main conclusions of this work are: 

1. All solid planets approximately follow the scaled 
mass-radius relationship log 10 R s — k\ + l/31og 10 (M s ) — 
k2M^ 3 for up to M p ~ 4M S . This relationship can be 
scaled to physical units by the values mi and r% given 
in Table 4. The corresponding planet mass, in physical 
units, to which the above mass-radius equation is valid 
ranges from M = 17M® to M p = 4OM0, depending on 
the material (see the mi values in Table 4.) 

2. There is no simple power law for the mass-radius re- 
lation; the mass-radius curve slope changes even within a 
relatively narrow mass range. 

3. We can use the same formalism described in the above 
equation to summarize the mass-radius relationships we 
computed for differentiated planets, using the scaling pa- 
rameters listed in Table 4. Given the uncertainties in the 
equations of state and the best expected future observa- 
tions, this simple, handy approximation supplies enough 
detail and accuracy to interpret any forecasted observa- 
tions. 

4. Highly detailed interior planet models are not needed 
to infer a solid exoplanet's bulk composition from its mass 
and radius. This is because the temperature structure and 
phase changes have little impact on the total planet mass 
and radius. 



a. Low-pressure phase changes (at < 10 GPa) are not 
important for a planet's radius because for plausible planet 
compositions most of the mass is at high pressure. For ex- 
ample, 97% of the planet's mass is at high pressures (>3 
GPa). For high pressure phase changes we expect the as- 
sociated correction to the EOS (and hence derived planet 
radii) to be small because at high pressure the importance 
of chemical bonding patterns to the EOS drops. 

b. Temperature can be approximated as a thermal pres- 
sure term. This thermal pressure causes a decrease in den- 
sity of on order three percent or less at relevant temper- 
atures and pressures. A change in average density of a 
planet translates into a smaller change in the planet ra- 
dius, because R p ~ Pp^ 1 ^ 3 ■ Conceptually, At low pres- 
sures (< 10 GPa) in the outer planetary layers, the crystal 
lattice structure dominates the material's density and the 
thermal vibration contribution to the density are small in 
comparison. At high pressures the thermal pressure con- 
tribution to the EOS is small because close-packed nature 
of the materials prevents structural changes from thermal 
pressure contributions. 

5. We identified several interesting properties of exo- 
planets. 

a. Planets are not likely to be found that have radii 
smaller than a pure Fe planet. Fe is the most dense el- 
ement out of which planets are expected to form. While 
not a new conclusion, this point is useful to keep in mind 
for designing and interpreting radius observations of exo- 
planets. 

b. Planets above the H2O curve must have a signifi- 
cant H/He envelope. We can therefore easily distinguish 
between exoplanets with significant H/He envelopes and 
those without, as is the case for GJ 436b. We therefore 
define a "super Earth" to be a solid planet with no signif- 
icant gas envelope, regardless of its mass. 

c. Because of their unique position on the mass-radius 
diagram, H2O planets with more than 25% water ice by 
mass can be identified with approximately 5% fractional 
uncertai nty in M p and R v . (A s imilar conclusion was also 
found bv l Valencia et alj (|2007bh ; that if the 7.5 M e planet 
they modeled has a large water content it would be identi- 
fied by a large radius.) Discovering a water planet orbiting 
at a small semi-major axis would serve as strong evidence 
for planet migration, since presumably water planets form 
more efficiently far from their host stars, beyond the ice 
line. This point is valid provided the planet is not a carbon 
planet and provided there is no substantial contribution to 
the planet radius from an atmosphere or envelope. 

d. Even with planet mass and radius measurement un- 
certainties better than 1%, planets of different interior 
composition can have the same total mass and radius. In 
other words, different mass fractions of iron cores, silicate 
mantles, and water outer layers can have the same total 
radius for a planet of the same mass. 

e. Carbon planets, if they exist, have mass-radius re- 
lationships that overlap with mass-radius relationships of 
non-carbon planets (i.e., water and silicate planets). This 
is because the zero-pressure density of graphite is simi- 
lar to that of H 2 and the zero-pressure density of SiC 
is similar to that of MgSiOa. Exoplanet atmospheres of 
transiting exoplanets will have to be observed and studied 
to distinguish between carbon planets and water/silicate 
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planets. 

6. We conclude that, while detailed interior structure 
models are needed to understand the atmosphere forma- 
tion and evolution, detailed interior structure models are 
not needed to infer bulk composition from exoplanet mass 
and radius measurements. 

We dedicate this study of mass-radius relationships 
to the memory of our co-author Cathy Hier-Majumder. 
It is with great sadness that we acknowledge her un- 
timely death. We thank Rus Hemley, Y. Fei, and Dan 



Shim for extremely useful discussions about high pressure 
physics and equations of state. We thank Mercedes Lopez- 
Moralez, Dimitar Sasselov, Jim Elliot for useful discussions 
about observational uncertainties. We thank Diana Va- 
le ncia for providing a m ore detailed version of Figure 4 
in IValencia et alj (|2006t ). as well as for interesting discus- 
sions about unpublished work. This work was supported 
by the Carnegie Institution of Washington, the NASA As- 
trobiology Institute, and the Massachusetts Institute of 
Technology. 
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Table 1 

Parameters for th e Vin et (V) or Birch-Murnagh am ( BME) EOS fits. Refer ence s: (T) IAhrensI (12000T). (2 ) 



IHanfland et al.I (11989T ); ( 3) ITakahashi fc SpainI (119891 ) ; (4) IAnderson et al.I (l2001fl: (5) 



IHemley ET AL.I (1198711: f7lHALLipAY ET AL.I (120031) 
(|2000ft ; (11) IKnittle 



. King 

(8) IDuffy et al I (11995F ): (9) IOlingerI (|197^ 



PrewittI (119891) (6 ) 



_ , , .. . ,, (10) IKarki et AL.I 

JeanlozI (jl987l ): (12) IAlexsandrov et al.I (jl989l ). a) seawater at 12°K. b) K' values are 

ASSUMED. C) A FOURTH ORDER BME FIT WAS USED WITH Kq = -0.016/GPA. 



V (cm^/mol) 


P (kg m 6 ) 


P (GPa) 


10.998300 


1.636617 


2.320 


10.429585 


1.725860 


4.155 


9.880818 


1.821712 


6.664 


9.351623 


1.924800 


9.823 


8.350443 


2.155574 


18.791 


7.878082 


2.284820 


25.361 


7.423411 


2.424761 


33.744 


6.986806 


2.576285 


44.314 


6.567891 


2.740606 


56.970 


6.165913 


2.919275 


74.188 


5.780497 


3.113919 


94.406 


5.411641 


3.326163 


126.815 


5.164734 


3.485175 


155.924 


4.654734 


3.867031 


240.696 


4.195170 


4.290649 


351.114 


3.780772 


4.760933 


498.660 


3.407399 


5.282621 


691.938 


3.070912 


5.861450 


937.585 


2.767547 


6.503954 


1260.182 


2.494293 


7.216474 


1673.049 


2.248138 


8.006625 


2188.301 


2.026072 


8.884186 


2853.712 


1.825836 


9.858497 


3691.387 


1.645548 


10.938603 


4737.211 


1.483327 


12.134882 


6040.611 


1.336538 


13.467635 


7686.171 



Table 2 

Density functional theory (DFT) EOS for water ice VIII and X. DFT predicts a gradual transition between 

THE TWO PHASES. 
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Material po [kg m a ] c [kg m 3 Pa ™] n 



Fe(a) 


8300.00 


0.00349 


0.528 


MgSi0 3 (perovskite) 
(Mg,Fe)Si0 3 


4100.00 


0.00161 


0.541 


4260.00 


0.00127 


0.549 


H 2 


1460.00 


0.00311 


0.513 


C (graphite) 


2250.00 


0.00350 


0.514 


SiC 


3220.00 


0.00172 


0.537 



Table 3 

Fits to the merged Vinet/BME and TFD EOS of the form p{P) = po + cP n . These fits are valid for the pressure 

range P < 10 16 Pa. 



Material 


m x [Me] 


n 


-Re J 


Pi [GPaj 


fci 


k 2 


k 3 


Fe(a) [modified-polytropic EOSJ 

MgSi03 (perovskite) [modified-polytropic EOS] 

H 2 (ice) [modified-polytropic EOS] 


5.80 

10.55 

5.52 


2.52 
3.90 
4.43 


1192 

693 

114 


-0.209490 
-0.209594 
-0.209396 


0.0804 
0.0799 
0.0807 


0.394 
0.413 
0.375 



Fc(a) 


4.34 


2.23 


MgSiOs (perovskite) 


7.38 


3.58 


H 2 (ice) 


8.16 


4.73 


Fe(0.675)/MgSiO 3 (0.325) 


6.41 


3.19 


Fe(0.3)/MgSiO 3 (0.7) 


6.41 


2.84 


Fe(0.225)/MgSiO 3 (0.525)/H 2 O(0.25) 


6.41 


3.63 


Fe(0.065)/MgSiO 3 (0.485)/H 2 O(0.45) 


6.88 


4.02 


Fe(0.03)/MgSiO 3 (0.22)/H 2 O(0.75) 


7.63 


4.42 



Table 4 

Conversion factors for the scaling relationships for equation (|23|) . The conversion factors give the physical 
values from the scaled parameters mass, radius, pressure, and density. see equations (ft3|) through (fttjj) . the 
first three rows of this table additionally give parameters for mass-radius relationships computed from the 
modified-polytropic-eos; these include the fej values for a fit to equation (123]) valid for m a > 4. 
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Fig. 1. — Equations of state for different materials at zero or 300 K. The black solid line is the EOS used in this study. The green dotted 
line is a fit to experimental data (either Vinet (V) or BME (BM)), appropriate for low pressures (typically below 200 GPa). The red dashed 
line is the TFD EOS, appropriate for high pressures (typically above 10 4 GPa). We adopt the Vinet or BME EOS at low pressures and switch 
to the TFD EOS at high pressures. Note that the abrupt increase in density of the MgSi03 BME curve at 2.8 X 10 13 Pa is above the pressure 
where we switch over to the TFD EOS, and illustrates the invalidity of extrapolating the BME to high pressures. 
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4 8 12 16 

log 10 P (Pa) 

Fig. 2. — Equations of state for five different materials used in this study. Each EOS data set was derived by combining a fit to experimental 
results and the TFD limit at high pressure. For water ice we used density functional theory as a bridge between experiment and the TFD 
theory. The EOSs are all reasonably well approximated by a polytropic-like expression p(P) = po + cP n , where po i s the zero-pressure density 
and c and n are constants. Table 2 lists these constants for some materials. The cyan triangles show one such fit for the H2O EOS. 
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Fig. 3. — The H2O EOS used in this study (solid cyan curve). We use the BME fit to the experimental data (Hcmlcv ct al. 1987, magenta 
triangles) up to P = 44.3 GPa. At this pressure we switch to the density functional EOS (squares). At P = 7686 GPa we switch to the TFD 
EOS (red dashed curve). 
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Fig. 4. — Mass-radius relationships for solid planets. The solid lines are homogeneous planets. From top to bottom the homogeneous 
planets are: hydrogen (cyan solid line); a hydrogen- helium mixture with 25% helium by mass (cyan dotted line); water ice (blue solid line); 
silicate (MgSi03 perovskite; red solid line); and iron (Fe (e); green solid line). The non-solid lines are differentiated planets. The red dashed 
line is for silicate planets with 32.5% by mass iron cores and 67.5% silicate mantles (similar to Earth) and the red dotted line is for silicate 
planets with 70% by mass iron core and 30% silicate mantles (similar to Mercury). The blue dashed line is for water planets with 75% water 
ice, a 22% silicate shell and a 3% iron core; the blue dot-dashed line is for water planets with 45% water ice, a 48.5% silicate shell and a 
6.5% iron core (similar to Ganymede); the blue dotted line is for water planets with 25% water ice, a 52.5% silicate shell and a 22.5% iron 
core. The blue triangles are solar system planets: from left to right Mars, Venus, Earth, Uranus, Neptune, Saturn, and Jupiter. The magenta 
squares denote the transiting exoplanets, including HD 149026b at 8.14 ijfg and GJ 436b at 3.95 Note that electron degeneracy pressure 
becomes important at high mass, causing the planet radius to become constant and even decrease for increasing mass. 
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Fig. 5. — Mass-radius relationships for planets with R p < 4Rgj. The curves are as in Figure[4] blue are for water ice planets, red are for 
silicate planets, and green is for pure iron planets. Black error bars are shown for 2% and 10% uncertainty in planet mass and radius. Each 
panel shows a different mass range. The terrestrial-mass solar system planets are shown with blue triangles. The exoplanets Gl 876d and 
Gl 531c are shown with a magenta square in panel c; although the radii are not known the are shown to represent known low-mass exoplanets. 
The exoplanet GJ 436b is shown in panel d. 
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Fig. 6. — Interior structure of solid exoplanets. From top to bottom the curves in panels a, b, and c are for planets with M p = 50, 10, 5, 
and 1 Mgj respectively. Panel a: silicate planets with a 32.5% by mass Fe core and a 67.5% MgSiC>3 mantle. Panel b: as in panel a but for 
planets with a 70% Fe core and 30% silicate mantle. Panel c: interior structure for water planets with 6.5% Fe core, 48.5% MgSiOs shell, and 
45% outer water ice layer. Panel d: interior model for two different water exoplanets with the same planet mass and radius: M p = 6.0Mj and 
Ftp = 2.0M^. The soli d curve is for a mo del with layers in percentages by mass of Fe/MgSi03/H20 of: 17/33/50 (similar to the composition 
of the water planet in (Lcjjcr ct al. 2004)) and the dotted line for 6.5/48.5/45 (similar to the composition of Ganymede). 
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M/M e M/M e 

Fig. 7. — Mass-radius relationships of solid planets to illustrate phase effects. Panel a and b: from top to bottom the solid curves are 
for: pure MgSiOs perovskite planets; MgSiC>3 planets with a 32.5% by mass Fe(e) core; MgSiOs planets with a 70% Fe core; and pure Fe 
planets. The dotted lines just above or beneath each of the solid curves is the planet mass-radius relationship for the BME or Vinet EOS 
alone without the TFD EOS. The top dashed line shows the pure MgSiOs perovskite planet with a phase change to MgSiOs enstatite at 
pressures less than 10 GPa. The squares show "super Earths" composed of a 32.5 % FeS core by mass and a 67.5% mantle; the mantle itself 
composed of 90% (Mg,Fe)Si03 by mass mixed with 10% MgO. Panels c and d: from top to bottom the solid lines are water ice planets with 
layers in percentages by mass of Fe/MgSi0 3 /H 2 of: 0/0/100; 3/22/75; 6.5/48.5/45, 22.5/52.5/25. The dotted lines beneath each of the 
solid curves is the planet mass-radius relationship for the BME form of H2O ice, without using our adopted H2O EOS (but using our adopted 
EOSs for MgSiOs and Fe). The dashed line shows the pure water ice planet with liquid water below 10 GPa. The magenta triangles show the 
mass-radius relationship for water planets with 17/33/50, illustrating a degeneracy with the 6.5/48.5/45 water planet (see also Figure|6ji). 
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Fig. 8. — The effect of a H/He gas envelope on planet radius, showing that only a small mass of H/He gas contributes a large factor to the 
planet radius. The lines in this Figure are the same as those in Figure [4] The magenta squares are planets with a fixed core mass composed 
of a mixture of 30% Fe by mass and 70% MgSi03. The cyan squares show the core mass and radius only. Panel a shows planets with a fixed 
core mass of 5Mj and 10Mg) (from left to right). Panel b shows planets with a fixed core mass of 20M^, 50M^ and lOOM® (from left to 
right). These H/He planet radii are underestimates because they are for zero temperature; temperature would make the planets larger for a 
given mass (i.e, move the squares up and left in this figure). 
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log 10 (M/M e ) 

Fig. 9. — Mass-radius relationship for carbon planets and helium planets. The curves are the same as those in Figure[4] The carbon planet 
mass-radius relationships are shown for: carbon monoxide planets (green triangles); graphite planets with 30% Fe cores by mass and 70% 
graphite mantles (circles); and SiC planets with 30% Fe cores by mass and 70% SiC mantles (squares). Pure cold He planet mass-radius 
relationships are shown by open squares connected by a solid line. The solar system planets are shown by the blue triangles. 
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Fig. 10. — Interior structure of homogeneous planets of 0.1, 1, 10, and 100 Mj. The top row shows the fractional mass as a function of 
fractional radius. The middle row shows the density as a function of fractional radius. The bottom row shows the pressure as a function of 
fractional radius. 
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Fig. 11. — The pressure level that contains 97% of the planet radius or mass. Panel a: the pressure at which the planet is 97% of its total 
size, as a function of total planet radius. Panel b: The pressure at which the planet contains 97% of its total mass, as a function of planet 
mass. The blue dashed curve is for pure water ice planets, the red solid curve is for pure silicate planet, and the green dotted curve is for 
pure iron planets, and H2O ice VII. 
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Fig. 12. — Thermal pressure for Fe (s), (Mg,Fe)SiC>3, and H 2 0. See text for details. 
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Fig. 13. — Isotherms for MgSi03 (perovskite) . Note that the perovskite phase of silicate does not exist for the whole region of low pressures 
shown in panel b, which is therefore shown to illustrate the general properties of silicates. 
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Fig. 14. — Isotherms for Fe(e) and H2O ice VII. See text for details. 



32 



1.2 



• i— i 

CO 
'cd 

o 

CO 



1 - 



0.8 - 



0.6 - 



0.4 - 



0.2 - 



~\ — i i 1 1 1 1 1| 1 — i i 1 1 nij 1 — i i 1 1 nij 1 — i i 1 1 1 1 1| 1 — i i 1 1 1 1 1| 1 — i i i mi 

n = 0.513 _ ~ 

n = 0.528 

n = 0.541 










0.001 0.01 











I 











0.1 1 10 

Scaled Mass 



100 1000 



Fig. 15. — Dimensionless mass-radius relationships for different materials. The scaled mass and radius are given in equations (115 I t and )16t . 
These dimensionless mass-radius relationships depend only on the exponent, n, in the modified polytrope EOSs (n = 0.513 corresponds to 
the H2O modified polytrope, n = 0.528 corresponds to Fe, and n = 0.541 corresponds to MgSiOs). For M s < 4, the mass-radius relationship 
for all materials takes approximately the same functional form. The squares are the approximate analytical scaled mass-radius relationship 
given in equation H39I I. which is valid for M B <C 1 and is the same for all n. 
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Fig. 16. — Scaled interior structure models for solid exoplanets computed from the polytropic-like EOSs. Panel a: scaled mass vs. fractional 
scaled radius. Panel b: scaled density vs. fractional scaled radius. Panel c: scaled pressure vs. scaled fractional radius. In all panels the solid, 
dotted, and dashed curves are for scaled masses of 0.5, 1, and 5, respectively. Three different values of n are used 0.528 (Fe; black), 0.541 
(MgSiOs; red), and 0.514 (H2O; blue). In all panels the solid and dotted curves are shown for only a single value of n, because the curves for 
all n overlap. The scaled interior structure solutions deviate from each other for the M B = 5 model. 



